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Abstract 

We study theoretically the cooperative light emission from a system of N ^ 
1 classical oscillators confined within a volume with spatial scale, L, much 
smaller than the radiation wavelength, Aq = 2ttc/ujo. We assume that the 
oscillators frequencies are randomly distributed around a central frequency, 
uq, with some characteristic width, 0, <^ ujq. In the absence of disorder, 
that is r2 = 0, the cooperative emission spectrum is composed of a narrow 
subradiant peak superimposed on a wide superradiant band. When O / 0, 
we demonstrate that if A'^ is large enough, the subradiant peak is not simply 
broadened by the disorder but rather splits into a system of random narrow 
peaks. We estimate the spectral width of these peaks as a function of A^, L, 
Q, and Aq. We also estimate the amplitude of this mesoscopic structure in 
the emission spectrum. 
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I. INTRODUCTION 



The study of cooperative phenomena in optics was initiated by the pioneering work of 
Dicke.0 The underlying physics of the cooperative emission can be readily understood using 
a classical approach. Suppose, that a large number of identical oscillators with frequency 
ujq are confined within a small volume with characteristic size L ^ 27cc/ujq = Aq, where 
Aq is the radiation wavelengh; this is reffered to as a "point" sample.! If r is the radiative 
lifetime of an isolated oscillator, then according to Dicke, the eigenmodes of the system of 
oscillators consist of one mode with a short lifetime t/N <C r, and A^— 1 modes with lifetimes 
much longer than r [by a factor ~ (L/Aq)^]. Correspondingly, the emission spectrum of this 
system consists of superimposed broad (superradiant) and a narrow (subradiant) bands. 
The intensities ratio of these bands is determined by the details of the excitation. This type 
of lifetimes redistribution is caused by the interactions among the oscillators through their 
radiation fields. 

Certainly, the classical picture does not describe all aspects of the cooperative emission. 
In fact, the original work of Dickelll primarily addressed the time evolution of the radiation 
emission, provided that at the initial moment, t = 0, all the oscillators are coherently 
excited. For this situation, the classical picture helps in understanding that the radiation 
is released during a short time, t/N; understanding of the initial stages of the emission 
process (the delay time statistics) requires, however, a quantum description.SS The original 
treatment in Ref. |^ also ignored the dipole-dipole interactions, which give rise to a spread in 
the oscillators frequencies (dephasing). The question whether or not this dephasing would 
completely destroy the cooperative emission is very non-trivial and was addressed in a 
number of later works. 

In the previous considerations of cooperative emission, it was assumed that all os- 
cillators (atoms, molecules or excitons) have the same frequencies. Such a restriction was 
adequate for the experimental situation in both gases and single crystals. To the best of our 
knowledge, the only account of disorder in the frequencies of the oscillators was given in Ref. 
^, which addressed the transient behavior of the cooperative emission from molecular aggre- 
gates. The case of J-aggregates corresponds to a symmetrical arrangement of oscillators in 
a circle. The author^ treated the disorder within the perturbation theory and averaged the 
second-order correction to the decay rates (the first-order correction vanishes upon averag- 
ing) with a Gaussian distribution. The advantage of the work in Ref. ^ is that the nearest 
neighbors dipole-dipole interactions were taken into account exactly. The drawback is in 
the perturbative approach, which rules out certain qualitative physical effects (see below). 

Whereas Ref. H addressed a rather particular situation, the following general questions 
might be asked. Suppose that the oscillators frequencies are randomly distributed with a 
characteristic width Q. Obviously, as Q increases, it would eventually destroy the cooperative 
features in the emission spectrum. Then what is the critical magnitude of How does this 
magnitude depend on the parameters of the system A^, L, and Aq? What is the structure of 
the emission spectrum when disorder is smaller than critical? 

These questions have become not purely academic due to the recent advances in the field 
of laser-action in 7r-conjugated polymers.BcJ Some experiments provide a strong evidence 
for cooperative emission from an ensemble of excitons in these materials, for excitation 
intensities exceeding a certain characteristic threshold.EUll On the other hand, it is well 
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known that the films of 7r-conjugated polymers are strongly disorderecfl (in the absense of 
disorder, cooperative emission by polymer chain was considered in Ref. |l^). They contain 
impurities and defects which break the polymer chains into segments of relatively short 
conjugation length, with a distribution depending on the film qualityJll This has a direct 
effect on the exciton energy, Tiuj, since it has been found that Tiuj directly depends on the 
chain conjugation length.ll^ 

The questions formulated above are addressed in the present paper. We study here 
the effect of disorder on the cooperative emission spectrum of the system of classical os- 
cillators. We consider the situation of incoherent excitation, which is most relevant to the 
experiment. In contrast to Ref. ^, we are interested in the nonaveraged (but universall 
properties of the emission spectrum. In other words, our goal is to assess the mesoscopi^ 
aspects of the cooperative emission. By mesoscopic we mean that, in the presence of a 
disorder, the emission spectrum of a large number of oscillators develops a fine structure. 
The actual shape of this spectral structure represents the fingerprints of the distribution of 
the oscillator frequencies and positions for a given realization. At the same time, the char- 
acteristic period and amplitude of the fine structure are determined by the net parameters 
of the system: A^, L, and VL. 

The paper is organized as follows. In Section ^ we derive the expression for the emission 
spectrum of a system of classical oscillators coupled by their radiation fields. In Section 
m we study in detail a simplified model in which the coupling among the oscillators is 
independent of distance. The eigenmodes of a "point" sample in the presence of disorder 
are analyzed in Section The universal properties of the mesoscopic stucture in the 
emission spectrum for small and large (but still smaller than Aq) sizes L are discussed in 
Sections and 0, respectively. In Section [VI I| the effect of the dipole-dipole interactions 
is addressed. The conclusions are given in Section |V111. 



II. THE BASIC EQUATIONS 

We consider a system of A^ oscillators located at random points r^, with frequencies uji 
randomly distributed around a central frequency ujq with a characteristic width Vt. Each 
oscillator is driven by the radiation field E(r,t) produced by all oscillators. The equation of 
motion for the displacement Ui of a given oscillator i reads 

iij + = — rii ■ E(ri,t), (1) 
m 

where e and m are the dipole characteristics (effective charge and mass) and rij is a unit 
vector in the direction of the dipole moment. 

The current density, associated with the oscillators motion, can be written as 

J(r,t) = e^njMj5(r - Fj). (2) 

i 

The current J plays the role of a source, which generates the electric field E(r,t) according 
to 

1 •■ 47r • 
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where c is the speed of hght. 

Within the classical approach, the emission spectrum of the system should be calculated 
as follows. We assume that at the initial moment, t = 0, all oscillators are excited with 
different phases (pi, and that the radiation field at the initial moment is zero, E(r, 0) = 0. 
The evolution of E with time can be then obtained by solving Eqs. After taking the 

limit r — > oo and expanding the field into harmonics, the spectral intensity can be obtained 
as = |E(oo, ct;)p. 

To carry out this program, it is convenient to employ the Laplace transformation. The 
transformed functions Ui{p) and E(r, p) satisfy the following system of equations 

+ p^)ui{p) = —lii ■E{vi,p) + uo{pcos(j)i -Ui sin (4) 
m 

-.2 



AE(r,p) - — E(r,p) = —J^ni [p'^Ui{p) - Mo(pcos0i - cjjsin^j)] 5(r - r^), (5) 



where uq sin (pi and uJiUo cos (pi are the respective initial displacement and velocity of the 
ith oscillator. The solution of Eq. for E(r,p) can be presented as a superposition of 
eigenmodes, Ey(r), of the wave equation for the electromagnetic field, 

AE.(r) + 4E.(r) = 0, (6) 

where tu^ is the eigenfrequency. Assuming that the modes are normalized, / (irE^(r) = 1, 
we obtain the following expression for E(r,p) 

E(r,p) = -47re V \p'^Ui{p) - Uo{p cos (pi - Ui sin (pi)] ^^^-^^^^ E^(r). (7) 
Substituting Eq. (0) into Eq. (^), we get a system of coupled equations for the amplitudes 

/ 2 , 2n- / N 47re2 [n, ■ E^(ri)] [nj ■ E^(rj)] \ 2- f ^ ( , • , n 

(^i + V )Ui{p) = 2^ — p Uj{p) - Uo{p cos (pj - ujj sm (pj) 

+Uq {p COS (pi — uji sin (pi). (8) 

To simplify Eq. (|^), it is convenient to introduce new variables Vi{p): 

jP'Ui (p) 

Mp) = {p cos (pi - Ui sin (pi) . (9) 

Uq 

Then Eq. (|^) takes the form 

{u^ + p^)vi + SijVj = uji{ui sin (pi - p cos (pi), (10) 



where the coefficients 



^Anfl [n..E.(r.)][n,.E.(r,)] 
m ^ uj^+ p^ 
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describe the coupling between oscillators i and j via the radiation field. 

Let us now express the intensity, I{uj), in terms of the variables Vi{p). The expression 
for E(r,p) follows from Eqs. (|^ and {^) 

E(r,p) = -47reE^.(p)^^4^E.(r). (12) 

The Fourier transform of the electric field is obtained by replacing p by iuj in Eq. (|12]). In 
the limit r ^ oo, only the pole uj,^ = uj contributes to the sum over z/, so that 

E(r, u) 1,^00 [n. ■ K{ri)] E,{r)6{ul - u^). (13) 

iu 

This corresponds to taking the continuum limit for electromagnetic modes. The terms pro- 
portional to [E,^(r) ■ E^(r)], which appear after calculating |E(r, 0;)^ from Eq. (|T^), oscillate 
rapidly if 7^ z/. Therefore, only the terms with = u survive at large r. These terms 
contain products of the form [rij ■ Ei,(rj)] [n^ ■ Ei^{rj)]. Note that the same products enter 



into the coupling coefficients, Sij, defined by Eq. (|TT]). This allows us to present the final 



expression for the spectral intensity in a compact form 

I{uj)<xJ2v,{lmS^,)v;, (14) 

ij 

where Vi and Sij are calculated at p = iu. 

We assume that the spread of the oscillators frequencies due to the disorder is much 
smaller than the central frequency, Q <^ ujq. This means that the frequency dependence of 
the coupling constants is weak, so that Sijiiu) can be evaluated aX oj = ojq. The real part 
of Sij, which comes from the principle value of the sum over modes in Eq. (pT]), diverges 
for i = j. This divergency is the manifestation of the Lamb shift, well-known in quantum 
electrodynamics, and can be absorbed into Wj. At the same time, the imaginary part of Sa, 
which results from the pole = ujq, is finite. It determines the radiative lifetime, r, of an 
individual oscillator via the relation ImSaliuJo) = 2ujq/t. For a single oscillator in vacuum, 
the modes Ej, are simply plane waves, and the summation over u in Eq. ([TT|) recovers the 
textbook result 

r = (15) 

For i j, the coupling Sij between two oscillators depends on the ratio rj^/Ao, where rij 
is the distance between the oscillators, and Aq = 2nc/uJo is the radiation wavelength. For 
Tij ^ Ao, both real and imaginary parts of Sij oscillate rapidly with r^j, and the effect of 
coupling is negligibly small for a large ensemble of oscillators.El For Vij <^ Aq, the real part 
of Sij represents the dipole-dipole interaction of the oscillators i and j. It is convenient to 
present Sij in the form 

^ii = '^{Pij + ioiij), (16) 

where f5ij and aij are the dimensionless matrices of coupling between the oscillators, defined 

as 
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and 



a,; 



rij ■ n, 

^ 5 



An 



)(n, 



2r?- 



(18) 



Turning back to Eq. (|TOD, we note that since the distribution of oscillators frequencies is 
relatively narrow, that is f2 -C ujq, we can make some simplifications. Namely, for p = iu, 
the factor {ujf + p^) in the Ihs can be replaced by 2uJo{u!i — u), and the rhs can be written 
as —iLjQe~^'^\ Finally, after rescaling Vi by factor uJq, Eq. (p!0D takes the form 



(uJi - Uj)Vi + 



T 



-i<t>i 



(19) 



Equation (|T^ together with Eqs. (|T^) and (|T^)-(|TB|) allow us to calculate the spectral 



intensity I{uj) for any set of initial oscillators phases. For aij 



(z ^ j), the 



eigenfrequencies of the system are simply the frequencies of individual oscillators, and the 
the emission spectum represents a superposition of Lorentzian peaks centered at tUj. In the 
presence of nondiagonal coupling, the eigenfrequencies are those of cooperative eigenmodes 
which, in turn, are determined by the imaginary part of the coupling, aij. In the experiment, 
the measured spectrum represents the result of averaging over many excitation pulses.lli In 
order to simulate the experimental situation, we will assume the phases (pi to be uncorrelated 
random numbers and average the result for the spectral intensity over all (pi. 



III. A SIMPLE MODEL 

In this section we consider a simplified situation, in which Eq. (|T9|) with random frequen- 
cies UJi can be solved exactly and the expression for the spectral intensity can be obtained 
in a closed form. Following Dicke,0 we disregard the dipole-dipole interactions by setting 
(3ij = 0. Although this approximation is rather common, later on we will discuss it in more 
detail. Turning to aij, we note that since L'^/Xq -C 1, the second term in Eq. (|1^) is a small 
correction to the first term. We therefore approximate the non-diagonal elements of aij by 



replacing r^JX^ with its average. 



a 



an, ■ n 



aj; 



(20) 



where the coupling constant a, with a typical value (1 — a) ~ L'^/X^ <^ 1, is the same for 
all pairs. Note however, that the disorder coming from random orientations of n, is still 
included. Later we will use this model for the analysis of the system (|19D with realistic aij. 



A. General solution 



For the model coupling (20), the system of equations (|^) takes the form 
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Ui — u -\ — (1 — a ] 

T 



with vector s defined as 



Vi H — arij ■ s = — -e" 



N 



-i4>i 



S = 

1=1 



(21) 



(22) 



A closed equation for s can be obtained by multiplying Vi, found from Eq. (|2ll) , by rij and 
taking the sum over i. This yields 



n,; n,; ■ s 



T ^ uJi - uo + i{l - a)/T 2 ^ iOi - 00 + i{l - q)/t 



-i<t>i 



(23) 



Solving Eqs. ( pT]) and ( p3D for t>j and substituting the result into Eq. (0), we obtain for the 
spectral intensity after some algebra 



liu) oc — Im 



T 



-1 



where we introduced a function 



a vector 



and a tensor 



^M^(^) = E 



= E 



1 


Ui — uj + i{l — 


a)/T 






uji — CO + i{l - 


a)/T 







(24) 



(25) 



cijj — + i(l — a)/r' 



(26) 



(27) 



where nj„ are the components of rij 



B. Identical oscillators 



Let us first consider the case of identical oscillators having the same frequencies 
uoi = uoq, and dipole momenta all aligned in the same direction. Then we find from Eq. 



ooi — uj + i{l — a) / t' 
and after averaging over the initial phases 0j, we obtain from Eq. 

{N -l){l-a)/T , {l-a + aN)/T 



(28) 



I{uo) oc 



+ 



(cJo-cj)2 + (l-a)2/r2 (,^o_cj)2 + (l-a + aA^)2/r2 



(29) 
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The emission spectrum is a superposition of a wide and a narrow Lorentzians with spectral 
widths r ~ N/t and 7 = (1 — a)/r, respectively. In accordance to the classical result,0 the 
eigenmodes of the system of identical oscillators coupled via their radiation field represent 
a single superradiant mode with short radiation time t/N , and — 1 subradiant modes 
with radiation time much longer than that for an isolated oscillator, r/(l — a) ^ t. The 
superradiant mode is a symmetric superposition of oscillator states and is strongly coupled 
to the radiation field, whereas the coupling of the subradiant modes to the radiation field 
is suppressed. In this case, the frequencies of all — 1 subradiant modes are degenerate, 
and the spectrum consists of a single narrow peak of width 7 on top of much broader band 
of width r, as shown in Fig. 1. As can be seen, with decreasing 1 — a, the hight of the 
subradiant peak increases, whereas the amplitude of the superradiand band diminishes. 

C. Random frequencies 

Consider now the case when the oscillators frequencies are random, but orientational 
disorder is still absent, i.e. all dipoles are aligned in one direction. Again we have F^yiuj) = 
5fj,uf{uj), with f{ijj) = f'{uj) + if"{uj) given by Eq. ([25|). Then a straightforward evaluation 
of Eq. (p3) yields (after averaging over the phases) 



where the function fiiu) = f[{uj) + ifiioj) is defined as 



f/{(i-^r) + (f)'/fr 



-f" 



(30) 



/iH = Et VTT^- (31) 

In order to clarify the underlying physics, it is useful to express the spectral intensity in 
terms of the system eigenmodes. The eigenfrequencies u!k are determined by the equation: 

?Vv 

1 + -f{uu) = 0. (32) 
r 

Then the intensity Eq. (|30| ) can be simply rewritten as 

^M-i: ,.-^,.^ . (33) 

where d;^ = RetDfc is the eigenmode frequency and u'l = ImtD^ characterizes its width. Note 
that for Ui = uq, we have A^ — 1 degenerate eigenmodes with oj[ = uq, and Eq. (|33D turns 
into Eq. m). 



D. Disorder in orientations 

In the presence of the orientational disorder, the spectral intensity ( P^ depends, in 
principle, on the direction of each rij. However, for large A^, one can replace the product 
ni^riiy in Eq. (^Tf ) for F^^, with its average. 
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(34) 



Thus, we have Ff^^{uj) = |5^jy/(u;), so that the expression for the spectral intensity is similar 
to Eq. ( ^OD with the only difference that in the first term, the functions f{uj) and fi{uj) are 
now multiplied by 1/3. This results in a shrinkage of the superradiant emission band by the 
same factor. At the same time, the width of subradiant peak increases by a factor of 3. Thus, 
the orientational disorder has no qualitative effect on the cooperative emission spectrum. 
The reason is that the coupling (20) is separable, that is it depends on orientations via the 
product rij ■ nj. Furthermore, for realistic aij given by Eq. ([ISD , the main (first) term has 
the same separable form; therefore, the orientational disorder does not qualitatively affect 
the cooperative emission spectrum and will be disregarded in the rest of the paper. 



E. Numerical results 

In Fig. 2 we plot the normalized spectral intensity in the absence of coupling, i.e. a = 0, 
with increasing number of oscillators. Each spectrum is calculated with a computer gener- 
ated set of N random frequencies uji, which we have chosen, for simplicity, to be uniformly 
distributed in the interval (c^o — ^,uJo + Q). For convenience, the spectra corresponding to 
different are normalized and shifted in the vertical direction. It can be seen that the 
peaks are resolved in the spectrum as long as the disorder, Q, is larger than N/t. We also 
see that for sufficiently large N, the intensity peaks are washed out from the spectrum. 

In Figs. 3-6 we present the results for I{uj) calculated using Eq. ( ^0]) for several values 
of a close to 1. The striking feature of the emission spectrum is its mesoscopic character. In 
the presence of disorder, the narrow subradiant peak of Eq. (|29|) (see Fig. 1) is not smeared 
out due to a large spread in Ui, as in the case of uncoupled oscillators (see Fig. 2), but rather 
splits into a multitude of narrow peaks corresponding to the eigenmodes of the disordered 
system. Furthermore, although the curves are calculated with different random sets of 
frequencies, the overall pattern of the emission spectrum exhibits certain universal features. 
In particular, it can be seen by comparing Figs. 3-6 that with increasing A^, the random 
structure of the spectrum undergoes several transformations, and that the characteristic A^, 
at which the changes in the pattern occur, is sensitive to the proximity of a to 1. This 
indicates a rather non-trivial structure of the eigenmodes, which we address in the next 
section. 



IV. STRUCTURE OF EIGENMODES 

The eigenmodes of a system of A^ oscillators coupled through their radiation field are 
determined by the homogeneous part of Eq. (ITUI) (we set (3ij = in this section) 



{ui - uj)vi + - 51 = 0- (35) 



Since the typical values of (1 — aij) ~ rj^/Xl are small, we split the second term in Eq. ( |35[ ) 
into a sum of the main contribution, with = 1, and a correction proportional to {aij — 1). 
Analogously to the consideration in the previous section, we rewrite Eq. (|35D as 
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1 

{uji - uj)vi + + (Ti) = 0, (36) 
r 



with 

1 



Expressing Vi from Eq. (^) and taking the sum over i, we obtain 



i 1 + (Ti , . 

1 + -E ^ = 0- 38 

T ^ UOj — UO 
3 J 

The equation for o"i follows from substituting of fj, found from Eq. (p6D , into the definition 
of a„ Eq. (13), 

a„ + ^E^^(l + ^.) = 0. (39) 

3 3 



The solutions of Eqs. (^) and (|39D determine the complex frequencies of the eigenmodes, 



A. "Point" sample 

Let us first analyze the effect of disorder on a system with all aij = 1, corresponding to 
the limit of a "point" sample, i.e. (L/Aq)^ ^ 1. With cxj = 0, the real and imaginary parts 
of Eq. dD read 

-Y ^^l^ = (40) 

1 ^ LJ" 



T 



3 



{iOj - iO'y + iO'' 



This system of equations has two different solutions with a crossover between them governed 
by the parameter Qt/N. For large disorder, Q ^ N/t, it can be readily seen that only one 
term in each of Eqs. (^) and (^T|) contributes to the sum. In this case, the solutions are 



simply u = LUj + i/r, as if the oscillators were uncoupled. In fact, this conclusion could be 
anticipated. The above parameter represents the ratio of the mean frequency spacing (MPS) 
of oscillators, Q/N, and the inverse lifetime of an individual oscillator, 1/r; when the former 
is much larger than the latter, Q/N ^ 1/r, the oscillators do not "feel" each other. 

In the opposite case of large (or weak disorder), N ^ Qr, the analysis of Eqs. ( ^OD 



and (|41| ) is carried out as follows. First note that in Eq. (|40|), which determines the real parts 
of the eigenfrequencies, u'j^, all the terms in the sum contribute now. Let us drop uj"'^ in the 
denominator of Eq. ( ^0|) (this step will be justified below). Then we obtain that the solutions 
a)^ are given by the extrema of the polynomial P{uj) = — uj). These determine the 

frequencies of the — 1 suhradiant modes. At the same time, in Eq. (^T]), which determines 
the imaginary parts of the eigenfrequencies, uj'l, all the terms in the sum are positive, so 



10 



that one should keep only the term with Uj closest to cu^. Since (cu^ — ^j) ~ ^/N for this 
term, we obtain the following estimate for the width of the subradiant mode: u'^ ~ 7, where 

7 ~ T^f/N'^. (42) 

It can be seen that 7 is much smaller than the MFS (by the factor Qt/N ^1). This justifies 
neglecting u"'^ in the denominators of Eqs. (^Dj) and (^T]). 

The superradiant solution of Eqs. (|40|) and (^) corresponds to the case u" ^ Q. Then we 
readily obtain Co' = N"^ J2j and cD" = T ~ N/t. We see that, indeed, T/Vtr^ N/Qt > 1. 
Therefore, the superradiant band in the spectral intensity is not affected by the disorder. 

We therefore conclude that cooperative emission is not destroyed by disorder. The spec- 
trum of the system consists of a single superradiant and — 1 subradiant eigenmodes. For 
large N/flr, the subradiant modes are well defined, since their spectral widths are much 
smaller than the MFS. 



B. Limit of weak disorder 

In this subsection, we address a nontrivial question about the fate of cooperative eigen- 
modes when the disorder in frequencies vanishes. In this limit, Qt/N 0, all oscillator 
frequencies become equal, i.e. uji — > ujq. In the absence of cooperative couphng, aij = 
(i 7^ j), the eigenfrequencies of the system are those of individual oscillators with the en- 
ergy width much larger than the MFS, 1/r ^ Q/N, so that the spectrum of the system is 
degenerate. 

However, the situation is more complicated in the presence of cooperative coupling, 
aij 7^ 0. Consider the case of a "point sample", aij = 1. In this case, the width of 
subradiant modes is given by Eq. (^2]). Important is that although the MFS diminishes with 
decreasing Q, the width 7 decreases even faster: '-//{Q/N) ~ Qr/N —>■ 0. In orher words, 
in the presence of even a very weak disorder, the narrow subradiant peaks do not overlap. 
Therefore, the cooperative modes remain distinct even though the "bare" oscillator modes 
were already degenerate. 

In the case of general coupling, the width 7 of subradiant modes for small values of Qt/N 
will be determined by the fluctuations of a^, as we will see below. 



C. Fluctuations of Oij 

Let us turn to the case with realistic coupling a^j. The eigenfrequencies Cok should now 
be determined from Eq. (^), which in component form reads 

with Oiiuj) = o-'-{uj) + ic'lioj) satisfying Eq. (|3^) , or in component form 
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- 7 ^ - (c.,-c.O^ + ..--^ = °- 

For u" <^Q/N, the system (^3|)-(^) can be approximately solved in the same way as for a 
"point" sample. The corresponding condition will be derived in Section V. 

When evaluating the contribution to the Ihs of Eq. (0) coming from the first term in the 
numerator, one should keep only one term in the sum with Uj closest to u;^: {ujj —oj'^) ~ Q/N. 
Then we obtain 



"""''^ 'E-^h (47) 



UJi. ~ 



where we again dropped u""^ in the denominator. Since a- <^ a" ^ 1 (see Section V), the 
frequencies cj^ in Eq. (^) are the same as for the case aij = 1. Finding uf in the first order 
from Eq. (|^), and substituting the result into Eq. (PTj), we obtain 



1 — ^ C^ij 1 



(48) 



The second term is the sought correction to the width of the subradiant modes. Remarkably, 
this term turns to zero if the matrix elements atj are replaced by their average a. Indeed, 
in this case the double sum in Eq. (^) would factorize into a product of two sums, each 
vanishing due to the fact that cD^ are the solutions of Eq. (|^) (corresponding to aij = 1). 
Therefore, the widths of the subradiant modes are determined by the fluctuations, Saij, of 
the coupling parameters aij rather than the deviation of their average, a, from unity. It 
should be noted that this property is general: one can easily see by comparing Eqs. ( ^31) 
and (0) to Eqs. (|45|) and (|46|) that for aij = const, we have a" = and a- ^ 1, so that the 
eigenfrequencies Uk are unaffected. 



D. Discussion of the numerical results 

We are now in the position to explain the spectra shown in Figs. 3-6. For the model 
coupling: 1, fluctuations only in the diagonal elements are finite: 5aij = 

(1 — a)5ij. Substituting this 5aij into Eq. (|48|) [instead of (a^- — 1)] and keeping only the 
term with (cD^ — ujj) ~ fi/A^ in the remaining sum, we obtain 

(49) 

The above expression indicates that after the cooperative modes have been formed (at 
~ VLt), the system can be found in two different regimes characterized by the relative 
magnitude of the first and second terms in the rhs. For intermediate number of oscillators, 
VLt ^ N ^ rir {1 — the width decreaeses with increasing N, as can be seen by 
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comparing the bottom second and third curves in each of Figs. 3-6 (note that the lowest 
curves with N = 2 show no sign of cooperative emission). In this regime, the system behaves 
in the same way as a "point" sample. With increasing number of oscillators, the dependence 
on saturates, and the width is dominated by the fluctuations of aij. Correspondingly, the 
change in the pattern of the peaks in Figs. 3-6, calculated for different values of (1 — a), 
occurs at different A^, as can be seen by comparing the next two curves in each figure. Note, 
however, that with further increase in A^, the curves exhibit yet another change in pattern. 
Namely, the peaks get smeared out (the top two curves in each figure). This occurs when the 
value of (1 — a)/r exceeds the MFS, Q/N, which is inconsistent with the above analysis. The 
reason for such a discrepancy is that for large A^, the model with coupling independent of 
separation between oscillators, becomes inadequate, as we mentioned above. For the correct 
description of the peaks smearing at large A^, the spatial dependence of aij is crucial; this 
question is addressed in Section 0. Nevertheless, for A^ < f2r/(l — a), this model describes 
accurately the mesoscopic features of the spectral intensity, as shown in the next section. 



V. STRONG MESOSCOPICS REGIME 

Let us now estimate the typical width of the radiation eigenmodes due to the fluctuations 
in aij. Since the configurational average of the second term in Eq. (^ ) [with Saij instead 
of {aij — 1)] vanishes, we need to evaluate {{uj'lY). Using the fact that only diagonal terms 
in the average {6aij6ai>j/) survive and omitting the first term in Eq. (^8]), we write 

The sum is dominated by the terms with {ui — tD^) ~ {uj — cD^) ~ Q/N. Since the typical 
spatial separation between two oscillators with close frequencies is ~ L, the separation 



fluctuations are of the same order. Thus, the typical fluctuation of aij is 6a = ^ {{6aijY) 
[L/Xof ^ and we finally obtain the typical width of a subradiant mode, 7 = y {{oj'^y), as 



7~-~-(^-j. (51) 

Comparing Eq. (|5l|) to Eq. (|^), we see that fluctuations in aij dominate the width 7 for 
N > nriXo/L). 

In order to characterize the fine structure in the emission spectrum, it is convenient to 
introduce the dimensionless parameter 



n[ L 



(52) 



It represents the product of a small, Qt/N, and a large, {Xo/L)"^, factors, which characterize 
the disorder and the system size, respectively. In terms of k, the condition for the formation 
of the cooperative modes, N/VLt ^ 1, can be presented as k ^ (Xq/L)'^. 



Using Eq. (52), the width (^) can be expressed in terms of the MFS as 
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(53) 



This result applies when k < Xo/L. On the other hand, it was implicit in the above 
derivation (Section pV]) that typical and a" are smaller than unity . The latter parameters 
can be estimated in a similar way from Eqs. (|45| ) and ( ^61) with the result: a" ~ and 
cr' ~ 0""^ ~ K,~^. Thus, the lower boundary for k, at which Eq. (|53D applies, is k ^ 1. For 
1, all terms in Eqs. (il)-(i|) become of the same order of magnitude, and for smaller 
K this system has no subradiant solutions, as discussed above. 

Since the MPS exceeds the width 7 within the entire domain 1 < k < Xq/L, the fine 
structure in the spectral intensity I{uj) is well pronounced. In other words, this domain 
corresponds to the strong mesoscopics regime. The opposite case k < 1 is considered in the 
next section. 



VI. WEAK MESOSCOPICS REGIME 

In the domain k <^ 1, the system cannot sustain eigenmodes that involve all oscillators. 
As a result, the eigenmodes become localized, in the sense that each eigenmode would 
comprise some iVc -C oscillators and occupy the volume with characteristic size Lc <C L. 
The magnitude of and A'c can be estimated from the following argument. Let us divide 
the system of oscillators into subsystems of increasingly smaller size. When the size of the 
subsystem becomes ~ Lc, the system of equations ( ^31 ) and (0), applied to a subsystem, 
first acquires a solution. This happens when the width 7c ~ t^^{Lc/ Xo)^ , determined from 
Eq. (^Tj) for a subsystem, becomes of the order of MPS within a subsystem, i.e. 

- ~ -. (54) 

t\\J ^ ' 

Taking into account that = N{L(./ , we find 

Lc ~ K^/^L, Nc ~ K^/^N. (55) 
Substituting these results back into Eq. (Q), we find for the eigenmodes width 

Prom Eq. (pSf ), we can also estimate how the relative amplitude of mesoscopic fluctuations 
in the spectral intensity I{uj) falls off with decreasing k: 

It is apparent that smearing of the fine structure in the cooperative emission spectrum with 
decreasing k occurs rather slowly. 
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VII. DIPOLE-DIPOLE INTERACTIONS 



In this section we study the effect of dipole-dipole interactions on the cooperative emis- 
sion from a disordered system. Note that for a "point" sample with L <C Aq, the typical 
magnitude of the dipole-dipole interaction between two oscillators is much larger than their 
superradiant coupling, Pij/cnj ~ {Xq/LY ^ 1. The structure of the eigenmodes in the 
absenceofsuperradiant coupling, given by Eq. (|l|) with aij = 0, was considered in several 
papers.c3cZl Renormalizat ion-group arguments of Ref. ^ (see also Ref. 0) suggest that all 
eigenmodes are delocalized. Numerical studie^i^i indicate a wide range of spatial scales 
in eigenmodes and, thus, seem to support this conclusion. In Ref. the role of general 
random-matrix perturbation in the spectrum of multilevel system was studied analytically; 
he ensemble-averaged renormalization of the spectrum of the system was derived which does 
not capture, however, the mesoscopic effects. Below we argue that finite disorder in combi- 
nation with superradiant coupling lead to a certain "resistance" of the system to large, but 
zero on average, dipole-dipole terms because of the formation of cooperative modes. 

In the absence of superradiant coupling {aij = 0) the dipole-dipole interactions lead to 
the shifts in the frequencies of individual oscillators. The resulting additional spread in cuj 
is, in general, much larger than the "bare" spread fl. This can be readily seen from the 
lowest-order correction to the frequency, Sui, which has the form 

1 B^- 

5^^ = -,J:^^^ (58) 

UJi — UJj 

(since = 0, the lowest-order correction to uJi is quadratic). The main contribution to 
the sum comes from pairs of oscillators located closely in space, with Vij ~ LN~^^^ (nearest 
neighbor interaction), so that 

/\ \3 

(59) 




Since the typical frequency difference for such pairs is ~ f2, we obtain 

'^'-&(f!- (60) 

On the other hand, the results obtained in the previous sections apply only if the additional 
disorder, caused by dipole-dipole interactions, does not affect the MFS. This requires the 
condition 5uJi ^ i7 to be met. Using Eq. (|3), this condition could be rewritten as Xq/ L <^ 
[VLt /Ny/'^ . Since the formation of cooperative modes occurs only if VLt/N <^ 1, one could 
draw the conclusion that neglecting the dipole-dipole interactions would be inconsistent 
with our basic assumption L -C Aq. 

The resolution of this apparent contradiction lies in the observation that, in the presence 
of superradiant coupling, i.e. aij ^ 0, the true eigenmodes of the system are cooperative 
modes comprised from a large number of oscillators. Therefore, the relevant condition should 
involve the shifts, (JcD^, of the eigenmodes frequencies, rather than 5uJi. In the first order, 
Suj'f. is given by an expression similar to the second term in the rhs of Eq. (|i8| ) [with j3ij 
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instead of {aij — 1)]. Since this term vanishes on average, as discussed above, the typical 



shift, 6uj' = J {{6uj'j^)'^), can be estimated from [compare with Eq. 




\2 



mn^i (E^-I^t— (61) 




There are two main contributions to the sum in the rhs. The first comes from the nearest- 
neighbor interaction with f3ij given by Eq. The second contribution originates from 
the pairs (ij) which are close in frequency; for such pairs, f3ij ~ {Xq/L)^. Both contributions 
turn out to be of the same order of magnitude, resulting in 

(62) 

This result is smaller than 6ui in Eq. (|60D by the factor N{N/Qr){Xo/ L)^ ^ 1. Such 
a dramatic difference illustrates the "resistance" of a coupled system of oscillators with 
disorder in frequencies to dipole-dipole interactions, as mentioned above. This property 
can also be qualitatively explained as follows. The dipole-dipole interaction between two 
subradiant modes can be viewed as an interaction between a mode and the electric field, 
E(r), created by the dipole moments of oscillators making up the other mode. Since the 
number of oscillators in a mode is large, their electric fields effectively cancel each other, so 
that the resulting net field, E(r), varies in space much slower than those of the individual 
oscillators. Note now that a slowly varying electric field couples only weakly to a subradiant 
mode. In fact, the suppression of the dipole-dipole interaction between subradiant modes 
has the same physical origin as their decoupling from the radiation field: had the electric 
field E been uniform, the cooperative modes would not interact at all with each other. This 
is the reason why the corrections to u'^. and cD^' vanish on average, and consequently the 
typical 60' and 7 are determined by the fluctuations of Sij. In contrast, the frequency shifts 
of individual oscillators are due to their interactions with the nearest neighbors, so that no 
cancellations occur. 

Thus we arrive at the condition Xq/L <^ (ilr)^^^, or, in terms of the parameter k, 

(63) 

This condition should be consistent with the condition for the formation of the cooperative 
modes, k -C (Xq/L)^. We see that both conditions are satisfied for sufficiently large N, i.e. 

^ (Xq/L)^. To account for different mesoscopics regimes, it is convenient to present Eq. 
(B^) in the form 

(64) 

Then n = 3, 4, and 5 correspond to the "point" sample, strong mesoscopics, and weak 
mesoscopics regimes, respectively. 
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VIII. CONCLUSIONS 



The main result of the present paper is that disorder in oscillators frequencies does not 
destroy the cooperative character of the emission from a "point" sample, as long as the 
MFS, Q/N, is smaller than the linewidth of an individual oscillator, r^^. In the opposite 
case, when Q/N ^ r^^, the spectrum represents a system of non-overlapping Lorentzians 
with the width r~^. 

It is convenient to characterize the disorder in terms of the dimensionless parameter 
K = QtXq/NL"^. Below we summarize our results for the characteristic width, 7, of the 
subradiant peaks (in units of Q/N) for different domains of k: 

7 = ^$(«:,L/Ao), (65) 
where the dimensionless function $ has the following asymptotes 

$ = '^(^^^ ' for ^ < K < ' s^^Pl^' 

$ = fi;^^, forl<K<— ^, strong mesoscopics, 

Ij 

$ = K^^/^, for K ^ 1, weak mesoscopics. (66) 

For K ^ {Xq/LY, the spectrum corresponds to uncoupled oscillators. 

In Section |T|, we presented an exact solution of a model with simplified (separable) 
coupling Eq. (^Ol). This model describes accurately the first two ("point" sample and strong 
mesoscopics) regimes in Eq. (^). It becomes, however, inadequate in the third (weak 
mesoscopics) regime, giving a instead of the correct k~^^^ dependence for the period of 
mesoscopic structure in the cooperative emission spectrum. 

Throughout the paper we have considered a three-dimensional system of oscillators. 
When the oscillators are confined to a plane, only the results for k ^ 1 should be modified. 
In this case, repeating the consideration of Section [V^, we obtain $ = Also for 

the relative magnitude of mesoscopic fluctuations, {6 1 /I), instead of Eq. (|57D we obtain 
(51/1) ~ K-^/\ 

Note finally that in experiments, such as photoexcited excitons in polymer films, the 
number of oscillators is governed by the excitation intensity.lii Thus, for a given disorder, 
the crossover from the strong mesoscopics regime {k > 1) to the weak mesoscopics regime 
{k < 1) can be simply achieved by increasing the excitation intensity levehil 
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FIGURES 

FIG. 1. Spectral intensity of N identical oscillators calculated from Eq. ( p9|) plotted vs. 
Alo = lo — uq foT N = 10, and a = (long-dashed line), a = 0.5 (dashed line), a = 0.8 (dotted 
line), and a = 0.9 (solid line). 

FIG. 2. Uncoupled oscillators: Spectral intensity I{uj) plotted vs. Aa; = a; — for several sets 
of random oscillator frequencies with Or = 5.0 and a = 0. 

FIG. 3. Coupled oscillators: Spectral intensity I{uj) calculated from Eq. ( ^0| ) for several sets of 
random oscillator frequencies with Or = 5.0 and a = 0.8. 

FIG. 4. Same as in Fig. 3, but for a = 0.85. 

FIG. 5. Same as in Fig. 3, but for a = 0.9. 

FIG. 6. Same as in Fig. 3, but for a = 0.95. 
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FIG. 3 
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FIG. 6 
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